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Abstract. The partonic transport model BAMPS (a Boltzmann approach to multiparton 
scatterings) is employed to investigate different aspects of heavy ion collisions within a common 
framework based on perturbative QCD. This report focuses on the joint investigation of the 
^ l' collective behavior of the created medium and the energy loss of high-pT gluons traversing this 

O , medium. To this end the elliptic flow and the nuclear modification factor of gluons in heavy ion 

collisions at 200 AGeV are simulated with BAMPS. 

The mechanism for the energy loss of high energy gluons within BAMPS is studied in detail. 
For this, purely elastic interactions are compared to radiative processes, gg ggg, that are 
^ ' implemented based on the matrix element by Gunion and Bertsch. The latter are found to be 

OO ' the dominant source of energy loss within the framework employed in this work. 

\o ' 
o 

I 1. Introduction 

■ Experiments at the Relativistic Heavy Ion Collider (RHIC) have established that jets with high 

transverse momenta are suppressed in Au + Au collisions with respect to a scaled p + p reference 
[U [2]. This quenching of jets is commonly attributed to energy loss on the partonic level as 
^ ■ the hard partons produced in initial interactions are bound to traverse the hot medium, the 

quark-gluon plasma (QGP), created in the early stages of the heavy ion collision. Another 
major discovery has been the strong collective flow of the created medium Oil], that is usually 
quantified in terms of the Fourier parameter V2 and in this context often referred to as elliptic 
flow. 

Due to the large momentum scales involved the energy loss of partonic jets can be treated 
in terms of perturbative QCD (pQCD) and most theoretical schemes attribute the main 
contribution to partonic energy loss to radiative processes O [U [3 El [9l [ini [H]- The bulk 
properties of the medium, such as the elliptic flow, are on the other hand usually investigated 
within hydrodynamical models. The comparison of hydrodynamic calculations to data indicates 
that the viscosity of the QGP is quite small [T^], possibly close to the conjectured lower bound 
^ = ^ from a correspondence between conformal field theory and string theory in an Anti-de- 
Sitter space |T5] . 

It is a major challenge to combine these two aspects, jet physics on the one and bulk evolution 
on the other hand, within a common framework. Recently the efforts to combine pQCD-based 
jet physics with hydrodynamic modeling of the medium have been intensified, for instance 
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results from hydrodynamical simulations are used as an input for the medium evolution in 
jet-quenching calculations (see [Hj for an overview) and as ingredients in Monte Carlo event 
generators [15]. However, these approaches still treat medium physics and jet physics in the 
QGP on very different grounds. Moreover, so far no schemes are available that cover the full 
dynamics of the interplay between jets and the medium. 

Partonic transport models might provide means to investigate bulk properties of the QGP 
and high-energy parton jets within a common physical framework automatically including the 
full dynamics of the system evolution. In previous publications jl6i [17] we have explored the 
capabilities of the transport model BAMPS (a Boltzmann approach to multiparton scatterings) 
with this goal in mind. 



2. The transport model BAMPS 

BAMPS [18[ [19] is a microscopic transport model aimed at simulating the QGP stage of heavy 
ion collisions via pQCD interactions consistently including parton creation and annihilation 
processes. Partons within BAMPS are treated as semi-classical and massless Boltzmann particles 
and at this stage the model is limited to gluonic degrees of freedom. Thus Nf = is understood 
throughout this work. The strong coupling is fixed to Og = 0.3. Please also note that the term 
"jet" as used throughout this paper refers to a single gluon with high energy that traverses the 
medium and does thus not fully coincide with the experimental notion. 

The interactions between partons are based on leading order pQCD matrix elements from 
which transition probabilities are computed. These are used to sample the interactions of 
particles in a stochastic manner |18] . The test particle method is introduced to reduce statistical 
fluctuations. For elastic interactions of gluons, gg -f-)- gg, we use the Debye screened cross section 
in small angle approximation '^^J':^^^ = — The Debye screening mass is computed from 

the local particle distribution / = f{p,x,t) via m|) = dcT^Us J (2-77^3 ^^cf, where do = 16 is the 
gluon degeneracy factor for = 3. 

Inelastic gg -f-T- ggg processes are treated via the Gunion-Bertsch matrix element [20] 

i2 _ 727r2a2s2 487ra^q^ 

(qi + ml)^ ki[(ki -q^)2+ ml] 



W^'-g9^999\ — 1^2 _|_™2\2 i-2 r/i. , „„2 l V^S ) \^) 



where qj_ and k^ denote the perpendicular components of the momentum transfer and of the 
radiated gluon momentum in the center of momentum (CM) frame of the colliding particles, 
respectively. Detailed balance between gluon multiplication and annihilation processes is ensured 

by the relation \M.gg^ggg\^ = dc \Mggg^gg\^ 



The Theta function in (1) is an effective implementation of the LPM (Landau, Pomeranchuk, 
Migdal) effect [21], that describes coherence effects in multiple bremsstrahlung processes. This 
interference effect cannot be incorporated directly into a semi-classical microscopic transport 
model such as BAMPS, so the cut-off (A^ — r) ensures that successive gg — )• ggg processes 
are independent of each other, r is the formation time of the gluon emitted with transverse 
momentum k± and Ag denotes the mean free path, i.e. the time between successive interactions, 
of the parent gluon. When comparing the formation time to the mean free path of the parent 
gluon special attention needs to be paid to the frames of reference. This ultimately renders the 
cut-off dependent on the boost j3 between the plasma rest frame and the center of momentum 



(CM) frame in which (1) is evaluated [T7] and numerically further complicates the calculations. 



3. Gluon jets in a static medium 

In order to obtain baseline calculations and to gain a better understanding of the energy loss 
mechanisms in BAMPS, we study the evolution of high energy gluons in a simplified brick setup. 
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i.e. the propagation within a static thermal medium with fixed temperature T. Also, for reasons 
of computing time, possible effects of the propagating jet on the medium are neglected in most 
computations presented in this section. For example, the mean energy loss per unit path length 
dE/dx is then calculated as follows (c = 1) 

where i denotes the interaction type {gg — )• gg, gg — )• ggg and ggg — )• gg) and i?* is the interaction 
rate for process i and (AE^) is the mean energy loss in a single collision of type i. 
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Figure 1. Differential energy loss of a gluon jet in a static and thermal medium of gluons at 
T = 400 MeV, T = 500 MeV and T = 600 MeV. 

Left panel: Elastic interactions only. Right panel: Including gg o ggg processes. 



The differential energy loss from elastic interactions as shown in the left panel of Fig. [T] 
exhibits the expected (see for an overview) logarithmic dependence on the jet energy 
E and the dominant quadratic dependence on the medium temperature T, ^[2^2 ^ 

C/jvra^T^ In ( ) , where Cn is the quadratic Casimir of the propagating jet, Cr = Ca = Nc 
for gluons. For T = 400 MeV and E = 50 GeV we find an elastic energy loss of ^\2^2 ~ 
1.2GeV/fm that increases to ^|2^2 ~ 2GeV/fm at E = 400 GeV. 

The energy loss caused by radiative gg — )• ggg interactions in BAMPS is much stronger and 
by far dominates the total differential energy loss (right panel of Fig. [1]) of a high energy gluon. 
The differential energy loss from gg — )• ggg is rising almost linearly with the energy, for example 
resulting in a total dE/dx w 32.6GeV/fm at E = 50 GeV and T = 400 MeV. 

The large differential energy loss in gg — )• ggg processes, however, is not governed by 
excessively strong cross sections. The individual cross sections increase only slowly with the 
jet energy as seen in the left panel of Fig. [2l For instance at E = 50 GeV and T = 400 MeV 
a gluon jet interacts with cross sections {(Tgg^gg) ~ 1.3 mb and {(^gg^ggg) ~ 3.5 mb. This 
emphasizes that BAMPS does indeed operate with reasonable partonic cross sections based on 
pQCD matrix elements. 

Rather it is the mean energy loss per single radiative process, (Aii^23)i shown in the right 
panel of Fig. [21 that determines the magnitude of the differential energy loss dE/dx. It is 
noteworthy that the strong and linear rise in the energy loss due to gg — )• ggg is only present 
when identifying the outgoing gluon with the highest energy as the jet gluon and thus using 
AE = E™ — max (S°"*, £^2"*) -^3"*) • This is the most natural choice and is employed throughout 
all calculations in this work. The average energy cj of the radiated gluon, however, is rising 
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Figure 2. Left panel: Cross sections for gg gg {2 ^ 2) and gg — )• ggg (2 — )• 3) processes. 
Right panel: Mean energy loss (Aii^23) in a single gg — )• ggg process. Upper group of lines: 
A^; = - max Lower group of lines (labeled "rad."): AE = u, with u 

being the energy of the radiated gluon. 

All quantities given as a function of the energy of the gluon jet inside a thermal medium with 
T = 400MeV. 

much slower with the jet energy. This is due to the fact that the energy is distributed among 
three outgoing particles, the gluon emitted with energy uj being only one of them. See [T7] for an 
in-depth discussion on how the complex interplay between kinematics, Gunion-Bertsch matrix 
element and phase space restriction imposed by the effective LPM cut-off gives rise to a fat tail 
in the AE23 distribution, that ultimately leads to large mean values {AE23). 

The full time evolution of the energy distribution of a jet particle propagating through 
the medium obviously contains more detailed information than the mean energy loss per 
unit path length. Fig. [3] shows p{E) dE\^, the probability that a parton starting with 
E{t = Ofm/c) = £'0 = 50 GeV has an energy E < E{t) < E + dE at a given time t. 




E [GeV] E [GeV] 

Figure 3. Time evolution of the energy distribution of a gluon jet that traverses a static and 
thermal medium of gluons (T = 400 MeV). The initial {t = Ofm/c) energy of the gluon jet is 
Eq = 50 GeV. 

Left panel: Elastic interactions only. Right panel: Including gg o ggg processes. 

For both gg — t- gg and gg — t- ggg the distribution of the jet energy induced by collisions with 
the constituents of the medium becomes rather broad. The distributions significantly differ from 
Gaussian shapes and a simple shift of the mean energy accompanied with momentum diffusion 
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could not account for the behavior depicted in Fig. [3l A distinct peak at lower energies only 
re-emerges at very late times. The mean energy loss as depicted in Fig. [T]is therefore a valuable 
observable but contains only limited information. It is noteworthy that there exists a finite 
probability for the jet to gain energy by collisions with the thermal gluons. 
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Figure 4. Left panel: Energy spectrum i^jj^^-^ of radiated gluons per energy interval du} and 
distance dx. to is the energy (lab frame) associated with the radiated gluon according to the 
Gunion-Bertsch matrix element. The energy of the gluon jet is = 50 GeV. 
Right panel: Angular distribution of the radiated gluon in the lab frame with respect to the 
original jet direction for different energies uj of the radiated gluon. Jet energy E = 50 GeV, 
medium T = 400 MeV. 



Accompanying the above discussions, Fig. [J] shows the energy spectrum of gluons radiated 
in gg — )• ggg processes and the angular distribution of the radiated gluons for different ranges 
of their energy uj. It is clearly visible that, due to the cut-off in transverse momentum (1) , the 
gluons cannot be emitted at very forward angles, an effect that is more pronounced for low uj. 
The spectra are peaked at energies uj <^ E, with a small tail reaching out to high energies. With 
increasing temperature the peak of the spectrum shifts towards higher energies in an apparently 
linear way, favoring the emission of gluons with higher energies. 



4. Au+Au collisions at 200 AGeV 

BAMPS has been applied to simulate elliptic flow and jet quenching at RHIC energies for 
the first time using a consistent and fully pQCD-based microscopic transport model to approach 
both key observables on the partonic level within a common setup. The left panel of Fig. 
shows that the medium simulated in the parton cascade BAMPS exhibits a sizable degree of 
elliptic flow in agreement with experimental findings at RHIC as established in [22\ 123] . And 
rj/s the gluon matter in BAMPS has been shown to be small |24j . 

For simulations of heavy ion collisions the initial gluon distributions are sampled according 
to a mini-jet model with a lower momentum cut-off po = 1.4 GeV and a X-factor of 2. The test 
particle method [18] is employed to ensure sufficient statistics and to allow for the resolution 
of adequate spatial length scales. The underlying nucleon-nucleon collisions follow a Glauber- 
model with a Wood-Saxon density profile and the results by Gliick, Reya and Vogt [25] are 
used as parton distribution functions. Quarks are discarded after sampling the initial parton 
distribution since currently a purely gluonic medium is considered. To model the freeze out of 
the simulated fireball, free streaming is applied to regions where the local energy density has 
dropped below a critical energy density Ec (ec = l.OGeV/fm^ unless otherwise noted). This 
setup has been successfully checked against experimental findings such as the distribution of 
transverse energy in rapidity and the flow parameter f 2 at various centralities in \19\ [22l [23] . 
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The right panel of Fig. [5] shows the gluonic Raa simulated in BAMPS for central, 6 = fm, 
collisions. It is roughly constant at R^^''^^ ^ 0.053 and in reasonable agreement with recent 
analytic results for the gluonic contribution to the nuclear modification factor Raa though 
the suppression of gluon jets in BAMPS appears to be slightly stronger. We expect improved 
agreement in future studies when employing a carefully averaged (6) that will be better suited 
for comparison to experimental data than the strict 6 = fm case. 
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Figure 5. Left panel: Elliptic flow V2 as a function of the number of participants for Au+Au 
at 200 AGeV for different combinations of the strong coupling and the critical energy density 
Ec- See [23] for more information. 

Right panel: Gluonic Raa at midrapidity (y e [—0.5, 0.5]) as extracted from simulations for 
central Au+Au collisions at 200 AGeV. For comparison the result from Wicks et al. [llj for the 
gluonic contribution to Raa and experimental results from PHENIX [26] for vr*^ and STAR [27] 
for charged hadrons are shown. 



As a first step towards making more extensive comparisons with experimental data and 
analytic models possible, we have computed the gluonic Raa for non-central Au + Au collisions 
at the RHIC energy of ^/s = 200AGeV with a fixed impact parameter b = 7fm, which roughly 
corresponds to 20% to 30% experimental centrality. The results is shown in the right panel of 

Fig. m 

A comparison in terms of the magnitude of the jet suppression for 6 = 7 fm is difficult since 
there are no published results for the gluonic contribution to Raa from analytic models available. 
Taking the ratio of the b = 7fm to the 6 = Ofm results as a rough guess indicates that the 
decrease in quenching is more pronounced in BAMPS compared to experimental data. The ratio 
of the nuclear modification factor between central (0% - 10%) and more peripheral (20% - 30%) 
collisions is Raa\o%-io% / Raa\2o%-30% ^ ^-^ experimental data, while for the BAMPS 

results RAA\b=Ofm/ ^AAlb =7fm ~ 0-4:. However, the issue of detailed quantitative comparison 
needs to be re-investigated once light quarks and a fragmentation scheme are included into the 
simulations. 

To complement the investigations of Raa at a non-zero impact parameter b = 7 GeV, we 
have computed the elliptic flow parameter V2 for gluons at the same impact parameter and 
extended the range in transverse momentum up to roughly pT ~ 10 GeV, see left panel of Fig. 
[6l For this calculation we have used a critical energy density Ec = 0.6GeV/fm'^ in order to be 
comparable to previous calculations. 

The V2 of high-pT gluons is rising up to ^ 4 GeV. Afterwards, from about pT ^ 5 GeV on, 
the elliptic flow slightly decreases again with pT- This behavior is in good qualitative agreement 
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Figure 6. Left panel: Elliptic flow V2 for gluons in simulated Au+Au collisions at 200 AGeV 
with b = 7fm. = O.GGeV/fm^. 

Right panel: Gluonic Raa as extracted from BAMPS simulations for b = Ofm and b = 7fm, 
Sc = l.OGeV/fm^. For comparison experimental results from PHENIX |26j for vr'' are shown for 
central (0% - 10%) and off-central (20% - 30%) collisions. 



with recent RHIC data [25] that for charged hadrons shows V2 to be rising up to V2 ~ 0.15 at 
~ 3 GeV followed by a slight decrease. 

5. Summary 

We have computed the gluonic contribution to the nuclear modification factor Raa and the 
elliptic flow V2 employing the pQCD based transport model BAMPS. This model provides means 
to investigate various characteristics of the evolution of the partonic medium created in heavy 
ion collisions, ranging from bulk properties to high-p^ physics, consistently including the full 
dynamics of the system. 

The gluonic Raa is found to be flat over a wide range in pT at Raa ^ 0.13 in off-central 
events (b = 7 GeV) and at Raa ~ 0.053 in central events (6 = OGeV) for a critical energy 
density of Ec = l.OGeV/fm^. Since BAMPS allows for the simultaneous investigation of high-p^ 
observables and bulk properties of the medium, we have also studied the elliptic flow parameter 
for gluons up to roughly 10 GeV for Au+Au at 6 = 7fm. V2 peaks at a ~ 4 + 5 GeV and 
slowly drops towards larger transverse momenta. 

In order to systematically investigate the energy loss of gluons as implemented in BAMPS 
we have studied the evolution of high energy gluons within thermal and static media of gluons. 
Inelastic gg — >■ ggg processes are found to be the dominant source of energy loss for high energy 
gluons in computations within the BAMPS model resulting in a strong differential energy loss 
that rises almost linearly with the jet energy. The strong mean energy loss in gg — t- ggg 
processes is due to a heavy tail in the AE distribution for single interactions, caused by the 
phase space configurations of outgoing particles dictated by the Gunion-Bertsch matrix element 
in combination with the effective LPM cutoff |17]. 

The characteristics of the strongly interacting, but still fully pQCD based, medium within 
the BAMPS description will be studied in further detail in upcoming works including light 
quark degrees of freedom. While a consistent modeling of low-pT hadronization needs 
careful consideration, the application of fragmentation functions to the high-py sector will be 
straightforward once light quarks are included and will allow for more direct comparison to 
hadronic observables. Also the application of BAMPS to heavy quark elliptic fiow and quenching 
will provide further valuable insight and is underway |29j . Additionally the medium response 
to high-pT particles will be studied in more detail. It has already been demonstrated (30j that 
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BAMPS offers the ability to describe collective shock phenomena in a viscous hydrodynamic 
medium. 
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